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The ocean rogue wave, one of the mysteries of 
nature, has not yet been understood or mod- 
elled satisfactorily, in spite of being in the in- 
tense lime-light in recent years and the concept 
spreading fast to other disciplines [l|-[sj- Rogue 
waves are extraordinarily high and steep surface 
waves. However, most of their theoretical models 
and exp erimental observations, excluding a few 
0. H. EEO- U] are one-dimensional, admitting lim- 
ited high intensity and steepness. We propose 
here a novel two-dimensional integrable nonlin- 
ear Schrodinger equation allowing an exact lump- 
soliton with special asymmetry and directional 
preference. The soliton can appear on surface 
waves making a hole just before surging up high, 
with adjustable height and steepness and disap- 
pear again followed by the hole. The dynamics, 
speed and the duration of the soliton is controlled 
by ocean currents. These desirable properties 
make our exact model promising for describing 
deep sea large rogue waves. 

RWs are reported to being observed in oceans, which 
apparently appear on a relatively calm sea from nowhere, 
make a sudden hole in the sea just before attaining sur- 
prisingly hi gh amplitude and disappear again without a 
trace jg, 1X2, |l3( . This elusive freak wave caught the imag- 
ination of the broad scientific community quite recently, 
triggering off an upsurge in theoretical 1(1 II , 14 , ll| and 
experimental [lH2( studies of this unique phenomenon. 
For identifying the class of such extreme waves the sug- 
gested signature of this rare event is a deviation of 
the probability distribution function (PDF) of the wave 
amplitude from its usual random Gaussian distribution 
(GD), by having a long-tail, indicating that the appear- 
ance of high intensity pulses are probable more often than 
that predicted by the GD [16]. In conformity with this 
definition RWs were detected in erbium doped fiber laser 
t| in optical lasers [l|, H[ , in nonlinear optical cavity 
i| , in acoustic turbulence in He II 0] and other set 
ups [l] . As for the theoretical studies on the ocean RW, 
apart from few nonlinear models based on the modula- 
tion instability, four- wave mixing and conformal mapping 
0, El, E3) or on some other effects [!, E(| , the nonlinear 
Schrodinger (NLS) equation 



admitting the Lax pair and exact soliton solutions [18| . 
Some models of RW generalise the NLS equation by 
adding extra terms on physical reasons like ocean current 



nonlinear dispersion [14|, [19j etc., which however 



id t q = d 2 2 q + 2\q\ 2 q, 



(1) 



is the most accepted one. Equation ([lj is a well known in- 
tegrable evolution equation in one space dimension (ID) 



makes the system non-integrable, allowing only numer- 
ical simulations. The most popular RW model used in 
many studies is a unique analytic solution of the original 
equation (fTJ, given by the Peregin breather (PB) or its 
trigonometric variant [1, 0, EH H, [HJ . The conventional 
soliton solution: q s = sechn(x — vt)e llykx+ult ^ of the NLS 
equation (p}, representing a localised translational wave 
behaves like a stable particle and unlike a RW propagates 
with unchanged shape and amplitude. In contrast, the 
exact PB solution of (p} 

q P {x,t) = e- 2tt {u + iv), u = G-l, v = -AtG, (2) 

where G — l/F, F(x,t) — x 2 + 4i 2 + j, represents a 
breather mode cos2i with unit intensity at both distant 
past and future, with its amplitude rising suddenly at 
t = to its maximum at x = 0, but subsiding again with 
time to the same breathing state. This intriguing RW 
like behaviour makes the PB a popular candidate for the 
rogue wave @, 0, El E | • 

Notice however that, the NLS equation ([1]) together 
with its different generalisations are equations in (1 + 1)- 
dimensions and therefore all of their solutions, including 
the PB, can describe the time evolution of a wave only 
along a line. Looking more closely into the PB we also 
realise that the maximum amplitude attained by the RW 
described by this solution is fixed, and just three times 
more than that of the background waves. The steep- 
ness of this wave as well as the fastness of its appearance 
are also fixed, since solution ([2]) admits no free param- 
eter. Therefore, though the well accepted PB or other 
solutions of the generalised NLS equation could be fit- 
ted into the working definition of the ocean RW, saying 
any wave with height more than twice the nearby signifi- 
cant height (average height among one-third of the high- 
est waves) could be treated as the RW [la ], they per- 
haps, with their severely restricted characteristics, can 
only explain moderately intense RW-like events in ID, 
as observed in water channels [3] , optical fibers 0, 0] or 
optical lasers [3, @] , but are grossly unfit as a model for 
the ocean RWs. First, ocean RWs are two-dimensional 
(2D) lumps appearing on the sea surface, which are not 
possible to describe by a ID equation like the NLS (pj. 
Second, the ocean RWs, as reported, might be as high 
as 17-30 meters in a relatively calm sea [1, E3], which is 
way above the value allowed by ID solutions like the PB. 
Third, the steepness of the oceanic RW and the speed of 
their appearance may vary from event to event, whereas 
they are fixed in the PB ([2]) admitting no free parameter, 
and are of values much below than those reported for the 
deep sea RWs. Note that in 2D water basin experiments 
as well as in the related simulations the amplitude and 
steepness of the RWs were found to be higher 1,@,E3, U| 



2 



than those predicted and observed in ID fA [Til] . 

The above arguments should be convincing enough to 
reject the PB together with other solutions of the gen- 
eralised ID NLS equation and look for a more realistic 
model for the ocean RWs. A straightforward 2D exten- 
sion of the NLS equation: 



id t q = did^q - d 2 dy V q 



2\q\\ 



(3) 



was proposed in some studies as a possible RW model 
l(l[llj]. However the 2D NLS equation © s not an inte- 
grable system and gives only approximate solutions with 
no stable soliton. Nevertheless, this unlikely candidate 
is found to exhibit RW like structures numerically, with 
higher intensity and steepness and with an intriguing di- 
rectional preference 0, [13] and broken special symmetry 

In the light of not so satisfactory present state in mod- 
elling the deep sea RWs, we propose a new integrable 
extension of the 2D NLS equation: 

id t q = d 1 dl 2 q-d 2 dy2q + 2iq(j x -j y ), j a = qd a q*-q*d a q 

(4) 

together with its exact lump-soliton as a suitable RW 
model. We find that, when the conventional amplitude- 
like nonlinear term in the non-integrable equation ([3]) 
is replaced by a current-like nonlinear term, the result- 
ing equation Q miraculously becomes a completely inte- 
grable system with all its characteristic properties, which 
is much rarer in 2D than in ID. More surprisingly equa- 
tion ((4]) admits an exact 2D generalisation of the PB with 
the desired properties of a real RW. Before proceeding 
further we notice that, our 2D NLS equation ((4]) can be 
simplified through a rotation on the plane to 



id t q + d 2 sv q + 2iq{qdyq* - q*d v q), 



(5) 



where the bar over the coordinates will be omitted in 
what follows. Encouragingly, we can find an exact stable 
soliton solution of the 2D equation ([5]) as 

q s (2d)(x,y,t) = sech«(y + px - vt ) e k ^+ k ^ t 

together with the associated Lax pair, infinite set of con- 
served charges and higher soliton solutions proving its in- 
tegrability (see supplementary Note) . However, the most 
important finding relevant to our present problem, is to 
discover an exact solution of ([5]) as a 2D extension of the 
PB. Before presenting the dynamical solution we con- 
sider first its static 2D lump-like structure localised in 
both space directions, describing a fully developed RW: 

qp(2d){x,y) = e 4lx (u + iv), u = G-l, v = -4xG, 

where G — — , F(x, y) = ay 2 + Ax 2 + c. (6) 
F 

One can check by direct insertion that (jSJ) is an exact 
static solution of the 2D nonlinear equation ([5]) with two 



arbitrary parameters a and c. It marks an important 
difference of solution ([6]) from the well known PB ([2]), 
which has no free parameter in spite of its close resem- 
blance with our solution. Looking closely into our 2D 
lump soliton ©, as shown in Fig. Id, we find that the 
wave attains its maximum amplitude: A rog = (l — c)/c — 
N — 1, at the centre x — 0, y — 0, where parameter 
c = 1/N is chosen through an arbitrary integer N. On 
the other hand, at large distances: |a;| — > oo, \y\ — > oo 
the wave goes to the background plane wave modulating 
as cos(4cc), with its amplitude decreasing to a constant: 
Aqo — 1, (Fig. la). The maximum amplitude relative 
to the background: ^p 2 - = N — 1, reachable by our RW 
solution, can be adjusted by choosing N to fit different 
observed heights of the RWs, and particularly the deep 
sea RWs having significantly higher amplitudes [t| [l3[ 
than usually assumed 

Eft 

Similarly, the steepness of 
the RW solution (j6]) at both sides as observed from the 
front: \d y qp^d) I, linked to another free parameter a, can 
also be changed to fit different observations. Note that 
the amplitude of the wave falls to its minimum: Aq = 0, 



at x = 0, y = ±yo with y = J i(l — c), which is rele- 
vant for the hole-wave formation as we see below. 

Our next aim is to construct a dynamical lump soliton 
out of the static solution ^ , to create a true picture of a 
RW which can appear and disappear fast with time. For 
constructing such a solution however we have to clarify 
first, whether it is possible in principle for our lump soli- 
ton to disappear, i.e. whether the soliton is free from all 
topological restrictions, which otherwise would prevent 
such a vanishing without a trace. The reason for such 
suspicion is due an interesting lesson from topology stat- 
ing that, when a complex field q{x, y) is defined on a 2D 
space with non- vanishing boundary condition |g| — > 1 at 
large distances, but having vanishing values q — > close 
to the centre, we can define a unit vector <b = A on an 

1-sphere S 1 . However, this vector tfi = (^1,^2) is well 
defined only at space boundaries: 9R 2 ~ S 1 (since q = 
at inner points), realising a smooth map: S 1 —¥ S 1 with 
possible nontrivial topological charge Q = n. This charge 
with integer values n = 0, 1, 2, . . ., labels the distinct ho- 
motopy classes and is defined as the degree of the map, 
which unlike a Nother charge is conserved irrespective 
of the dynamics of the system. Such a situation occurs 
for example in type II superconductors with the charge 
linked to the quantised flux of vortices for the magnetic 
field B(x,y) : 



2ttQ 



dS ■ B 



dl-A, 



(7) 



where B = curl A. = d x 4>id y <f>2 — d x 4>id y (f)\. Notice that, 
our complex field solution qp(2d){x 1 y) possesses clearly 
the features of <f> discussed above, since (j6]) goes to a con- 
stant modulation — e Alx at large distances and vanishes 
at points (0, ±yo)- Note that, such a solution related to 
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a sphere to sphere map can not go to a trivial configu- 
ration , if it belongs to a homotopy class with nontrivial 
topological charge: Q = N 7 N = 1,2,3,..., due to con- 
servation of the charge, with the only exception for the 
class with zero charge Q = 0. Therefore, for confirming 
the possible appearance/disappearance property of a RW 
for solution (|6]), we have to establish first that in spite 
of defining a nontrivial topological map, it belongs nev- 
ertheless to the sector with topological charge: Q = 0, 
i.e. our lump soliton is indeed shrinkable to the vacuum 
solution. For this we calculate explicitly the topological 
charge |7]) associated with ([6]) as 



2-kQ = J dl A = J (dxA x 



dyA x ),. 



(8) 



with A a = fadafa, 4>i = Req/\q\, 4> 2 = Imq/\q\, where 
the contour integral along x and y are taken along a 
closed square at the boundaries of the plane. Substitut- 
ing explicit form of solution q{x, y) from ([6]) and arguing 
about the oddness and evenness of the integrants with 
respect to x, y or checking directly by Mathematica one 
can show that the related charge is indeed Q = and 
therefore the solution belongs to the trivial topological 
sector as we wanted. The intriguing reason behind this 
fact is that, the two holes appearing here have in fact 
opposite charges resulting to the combined charge be- 
ing zero. For constructing a dynamical extension of the 
2D RW solution ^ we realise that, the sudden change 
of amplitude with time, as necessary to mimic the RW 
behaviour, would result to a non-conservation of energy. 
Therefore it can not be described by an integrable equa- 
tion alone, which insists on a strict conservation of all 
charges and therefore our integrable equation ((SJ needs 
certain modification. On the other hand, the importance 
of ocean currents in the formation of RWs, which is not 
considered in ((5]), is documented and repeatedly empha- 
sised [HI, 12 1 . This motivates us to solve both these 
problems in one go, by modifying our equation (|5|) with 
the inclusion of the effect of ocean currents, by adding a 
term in the form —ill c q v , as in [l3| . For obtaining an 
exact dynamical RW solution to our modified 2D NLS 
equation, we choose the currents flowing along transverse 
directions and changing with the time and location as 
U c (y, t) = As is apparent, the currents would flow to 
the centre of formation of the RW from both the trans- 
verse directions, with the speed increasing as they ap- 
proach, stopping however completely at the moment of 
the full surge, at t = 0. The picture gets reverted after 
the event with currents flowing back quickly, away from 
the centre. Though U c looks ill-defined, the multiplica- 
tive factor q y makes the current term well-behaved and 
the modified 2D NLS equation (([5]) with the inclusion 
of the current term) admits now an exact dynamical 2D 
Peregin soliton, which interestingly has the same form as 
(|6"|), only with the function F becoming dynamical with 



the inclusion of time variable as 

F(x,y,t) =F(x,y)+nt 2 , F(x, y) = Ax 2 + ay 2 + c. (9) 

The arbitrary parameter [i appearing in solution ([9]) is 
related to the ocean current and can control how fast 
the RW appears and how long it stays. It is convinc- 
ingly demonstrated in Fig. 1 , how this exact dynamical 
2D lump-soliton evolves from a background plane wave 
with a slight depression appearing already at the centre 
x = 0,y — (Fig. la). As the time passes a sudden 
hole is formed right in the centre (Fig. lb) at the mo- 

as told in marine-lore 
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ment t h = -Jj^O- ~ c), 
which splits into two and shift apart in the transverse 
direction (Fig. lc) to make space for a high steep up- 
surge of the lump forming the actual RW (Fig. Id) at 
time t = 0. The RW disappears fast into the background 
waves with the holes merging at the centre and vanishing 
again, thus describing vividly well the reported picture 
of the ocean RWs [H, H3, [l7j as well as those found in 
large scale experiments [3j. The surface waves modelled 
by our solution ([6]), as visible also from Fig. 1, show a 
distict directional preference and an asymmetry between 
the two space variables, 
a b 



□ 





Figure 1: Formation of 2D rogue wave and holes 

Panels a-c show the modulus of the dynamical 2D lump- 
solution and panel d its real part at time t = 0. Figures 
correspond to parameter values c = 1/6, a — 1.2, fi = 
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1.2. a, Almost constant amplitude background wave with 
a slight depression appearing at the centre (at time t = 
—30.00). b, Creation of a hole at the centre at time 
t h = —0.83. c, The hole splitted into two are drifting 
away from the centre (at t = —0.40). d, Full grown RW 
formed over the background modulation. The amplitude 
attained by the RW is five-times that of the background 
waves (linked to the choice of c). 

We conclude by listing a few distinguishing features 
of our proposed dynamical lump soliton (((6|) with 
which are important for a ocean RW model. 1) This is 
the first 2D dynamical RW model given in an analytic 
form. 2) It is an exact Peregin like lump soliton solution 
of a novel 1 + 2-dimensional nonlinear integrable equa- 
tion. 3) The dynamics of the solution is controlled by 
the ocean currents. 4) The height and steepness of the 
RW can be adjusted by independent free parameters. 5) 
The fastness of its appearance and the duration of stay 
can also be regulated by another parameter linked to the 
ocean current, which becomes the key factor in creating 
the dynamical RW. 6) The proposed solution exhibits 
a clear broken special symmetry in variables x and y as 
well as a directional preference, which are suspected to be 
crucial features in the formation of 2D RWs [!, 0, E3, El ■ 
6) Strange appearance (and disappearance) of a hole just 
before (and after) the formation of the rogue wave [IlT[l7j 
is also confirmed in our model. 

We hope that, this significant breakthrough in describ- 
ing large ocean RWs through a dynamical analytic lump- 
soliton would also be valuable for experimental findings 
of 2D RWs in other systems like capillary fluid waves |6[ 
optical cavity waves [J] and basin water waves Q . 
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Supplementary Note: Integrable properties 
of the Novel 2D NLS Equation (equation (5) in the 
main text) 

I. Lax pair representation: 

The integrable nonlinear equations can be linked to a 
linear system 

<S> X = $ t = V$, (10) 

associated with a pair of Lax martices U, V . The exis- 
tence of a Lax pair is considered to be a guarantee for the 
integrability of a nonlinear partial differential equation. 
The Lax pair dependent on an additional parameter A, 
called spectral parameter, are constructed in such a way 
that their flatness condition : 

u t - v x + [U, V] = 0, (11) 
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which is the compatibily condition of the overdetermined 
linear system, yields the given nonlinear equation 

For our nonlinear evolution equation (equation (5) in 
the main text) in 2D space we can find the Lax pair given 
through the standard 2x2 Pauli matrices a a , a = +, — , 3 
as 

U = t{U n <T 3 + U 12 a+ + U 2 i<r~), (12) 
where U\\ — 2A 2 + \q\ 2 , Uu = f/|i = 2Ag - iq y , and 
V = i(Vn<7 3 + V 12 <t+ + V 21 a-) (13) 

where 

V 11 = -4\ 3 + 2\\q\ 2 +i(qq: j -q*q y ), 

V12 = V 2 \ = -4A 2 9 + 2i\q y + \q\ 2 q + q yy , 

where the notation q y = d y q etc. as popular in the soli- 
ton commulity has been introduced. It can be checked 
directly that our integrable equation (5) (in the main 
text) is obtained from the flatness condition of the Lax 
pair mini). 

II. Conserved charges & soliton solutions : 

The infinite set of conserved charges associated with 
our integrable system can be derived systematically from 



the Riccati equation obtained from the first equation in 
(flC)]) with the Lax operator (IT2|) [2] . Exact higher soliton 
solutions of equation (5) (in the main text) can be ob- 
tained by the Hirota's bilinearization method, similar to 
the NLS solitons 



III. 2D rogue waves as exact solutions : 

The static 2D rogue wave (6) (in the main text) as 
an exact solution of the 2D NLS equation (5) (in the 
main text) can be checked by direct insertion. The The 
dynamical 2D rogue wave solution (equation (6) with (9) 
in the main text) can be checked as an exact solution, by 
inserting it in the modified 2D NLS equation obtained 
by adding the ocean current term —iU c q y to equation 
(5) (see the main text). 
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